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ABSTRACT 

The  reorder  point  and  reorder  quantity  for  a  multi- 
echelon  inventory  system  consisting  of  two  levels  were 
determined  through  the  use  of  a  mathematical  model  and  a 
computer  simulation  was  used  to  verify  the  results. 

The  main  echelon  supported  two  lower  echelon  stock 
points  and  reordered  using  a  continuous  review  inventory 
policy.   The  two  lower  echelon  stock  points  operated  under 
periodic  review  policies. 

The  measure  of  effectiveness  used  was  to  minimize  total 
system  costs  subject  to  a  constraint  on  the  maximum  number 
of  back  orders  per  year. 

The  results  from  the  mathematical  model  were  used  as 
input  values  for  the  simulation  and  measures  of  effective- 
ness were  compared.   An  alternative  procedure  was  proposed 
and  simulated;  and  the  results  of  the  three  products  were 
compared. 
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I.   INTRODUCTION 

Despite  the  fact  that  all  military  supply  systems  and 
many  large  civilian  corporations  have  multi-echelon  and 
multi-location  inventory  systems,  few  multi-echelon  inven- 
tory models  are  currently  being  used.   The  policies  that 
have  been  used  are  the  single  echelon  and  single  location 
inventory  policies  that  attempt  to  minimize  the  total 
variable  cost  of  a  single  location  ignoring  total  system 
cost.   The  difficulties  of  optimizing  a  multi-echelon  and 
multi-location  inventory  system  are  most  likely  due  to  the 
complexity  of  demand  and  the  interdependency  among  the  units 
in  different  echelons. 

A  mathematical  model  for  a  multi-echelon  and  multi- 
location  inventory  system  is  developed  in  Chapter  II.   The 
objective  of  this  model  is  to  minimize  the  total  annual 
variable  costs  subject  to  a  maximum  allowable  number  of 
back  orders  per  year  for  the  entire  system.   This  model 
consists  of  one  first  echelon  location  operating  under  a 
continuous  review  policy  and  two  second  echelon  systems 
which  operate  under  periodic  review  policies.   An  example 
is  presented  to  obtain  numerical  results. 

In  Chapter  III,  a  computer  model  is  developed  to  simu- 
late the  multi-echelon,  multi-location  system.  This  simu- 
lation model  was  used  to  check  the  results  obtained  from 
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the  mathematical  model  and  to  demonstrate  the  interaction 
between  echelons.   Graphical  plots  are  generated  which  show 
the  inventory  position  of  each  entity  in  the  model. 

In  Chapter  IV,  comparisons  are  made  of  the  simulation 
results  and  the  results  obtained  from  the  mathematical 
model.   The  comparisons  suggest  some  modifications  to  the 
operating  policy. 

The  system  was  simulated  again  using  the  modified 
operating  policy  and  substantial  improvements  in  the 
effectiveness  of  the  multi-echelon  system  were  observed. 
These  results  are  described  in  Chapter  V. 

Chapter  VI  summarizes  the  results  of  the  research  and 
concludes  with  some  suggestions  for  additional  research. 


II.   A  PROBABILISTIC  MATHEMATICAL  MODEL  FOR  MULTI-ECHELON 
INVENTORY  SYSTEM 

A.  DESCRIPTION  OF  MODEL 

In  this  model  there  are  three  systems.   One  of  these 
systems  represents  the  highest  echelon  and  is  called  the 
Main  System.   The  other  two,  at  separate  locations,  are 
the  lower  echelon  and  are  called  System  One  and  System  Two, 
respectively.   Each  system  carries  its  own  inventory,  and 
receives  random  demands.   System  One  and  System  Two  are 
dependent  upon  the  Main  System  but  independent  of  each 
other.   System  One  and  System  Two  can  only  be  resupplied 
by  the  Main  System.   In  other  words,  they  cannot  order 
from  any  other  suppliers.   The  Main  System  replenishes 
stocks  by  placing  orders  to  external  suppliers. 

B.  OBJECTIVE  OF  MODEL 

The  objective  of  this  model  is  to  minimize  the  total 
expected  annual  variable  costs  of  three  systems  subject 
to  a  specified  expected  number  of  total  back  orders  per 
year.   In  fact,  this  is  the  same  as  minimization  of  total 
yearly  variable  costs  subject  to  a  specified  minimum  level 
of  customer  satisfaction. 

Based  on  these  objectives,  decision  variables  for  each 
system  will  be  calculated  to  achieve  optimum  levels  for 
the  entire  multi-echelon  supply  system. 
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C.   ASSUMPTIONS  OF  THE  PROBABILISTIC  MATHEMATICAL  MODEL 
1.   Main  System 

It  is  assumed  that  the  Main  System  under  considera- 
tion consists  of.  a  single  installation  which  utilizes  trans- 
action reporting.   This  system  is  governed  by  a  (Q,  r)  type 
policy  with  back  orders. 

The  assumptions  in  addition  to  the  continuous  review 
assumption  are: 

a.  The  cost  of  operating  the  information  processing 
system  is  independent  of  Q  (reorder  quantity)  and  r  (reorder 
point) . 

b.  The  unit  cost  C  of  the  item  is  a  constant  in- 
dependent of  Q. 

c.  The  back-order  cost  is  constant  (II)  ,  per  unit 
back  ordered  regardless  of  the  length  of  time  the  back - 
orders  exist. 

d.  There  is  never  more  than  a  single  order  out- 
standing.  This  assumption  implies  that  when  the  reorder 
point  is  reached,  there  are  no  orders  outstanding;  there- 
fore, the  inventory  position  is  equal  to  the  net  inventory. 
Thus,  the  reorder  point  will  be  the  same  regardless  of 
whether  it  is  based  on  the  inventory  position  or  net 
inventory. 

e.  Procurement  lead  times  are  independent  and 
identically  distributed  (i.i.d)  random  variables  with  a 
gamma  distribution. 
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f.  All  variables  are  treated  as  continuous. 

g.  The  demands  are  Poisson  distributed  with  the 

mean  number  of  demands  per  year    a  constant  A  . 

^  m 

h.   The  reorder  point,  r,  which  is  based  on  the 
inventory  position  is  positive. 

With  a  back  orders  constraint,  it  is  infeasible  to 
wait  until  back  orders  exist  before  placing  an  order. 
Because  of  this  and  assumption  (d)  there  will  be  no  back 
orders  outstanding  at  the  reorder  point.   As  was  discussed 
in  assumption  (d) ,  at  the  reorder  point,  the  inventory 
position  is  equal  to  the  on-hand  inventory. 

For  this  model,  any  one  of  the  three  inventory 
levels  on  hand,  net,  or  inventory  position  can  be  used  to 
define  the  reorder  point;  and  the  reorder  point  has  the 
same  value  for  any  one  of  them.   It  should  also  be  noted 
that  to  use  the  on-hand  level  it  must  be  assumed  that 
after  an  order  arrives,  it  is  sufficient  to  fill  all  back 
orders  and  raise  the  on-hand  inventory  level  above  the 
reorder  point.   If  this  ever  failed  to  happen,  the  reorder 
point  would  never  be  reached  again  and  the  system  would 
continue  to  accumulate  back  orders  indefinitely. 

When  the  reorder  point  is  thought  of  in  terms  of 
the  inventory  position  of  the  system,  then  assumption  (d) 
guarantees  that  the  on-hand  inventory  will  always  exceed 
the  reorder  point  when  an  order  arrives;  otherwise,  it 
would  not  be  possible  to  have  only  a  single  order  outstanding. 
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2.   System  One 

This  system,  which  is  at  the  lower  echelon,  is 
governed  by  a  periodic  review  policy.   The  operating 
doctrine  is  the  most  widely  used  type  of  periodic  review 
which  is  the  order  up  to  R  policy.   All  demands  which 
occur  when  the  system  is  out  of  stock  are  back  ordered. 

For  this  periodic  review  system  the  time  between 
reviews  will  be  denoted  by  T,  and  at  each  review  time  a 
sufficient  quantity  is  ordered  to  bring  the  inventory 
position  of  the  system  up  to  a  level,  R,  regardless  of  the 
amount  of  on-hand  inventory.   This  policy  dictates  that 
at  review  times  even  if  the  inventory  position  is  R-l, 
only  one  item  must  be  ordered  to  bring  the  inventory 
position  up  to  R,  ignoring  the  high  cost  of  placing  the 
order.   Despite  its  appearance  of  being  illogical,  this 
assumption  simplifies  the  formulation  of  the  system  and 
is  very  unlikely  to  occur  on  high  demand  items.   When  an 
item  experiences  zero  demand  in  a  review  period,  there  is 
no  need  for  order  because  the  inventory  position  is 
already  at  R. 

The  other  assumptions  are: 

a.  The  cost,  J,  of  making  a  review  is  independent 
of  the  variables  R  and  T. 

b.  The  unit  cost,  C,  of  the  item  is  constant  and 
independent  of . the  quantity  ordered. 
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c.  Back  orders  are  incurred  only  in  very  small 
quantities.   This  implies  that  when  an  order  arrives,  it 
is  almost  always  sufficient  to  meet  any  outstanding  back 
orders. 

d.  The  backorder  cost  is  constant,  IT,  per  unit 
back  ordered  regardless  of  the  length  of  time  the  back 
order  exists. 

e.  Procurement  lead  times  are  i.i.d  random 
variables  with  a  gamma  distribution. 

f.  Orders  are  received  in  the  same  sequence  in 
which  they  were  placed.   It  should  be  noted  that  for  (Q,  r) 
models,  the  two  assumptions  that  orders  were  received  in 
the  sequence  placed  and  that  lead  times  are  i.i.d  random 
variables  could  not  both  hold  rigorously,  since  there 
exists  a  positive  probability  that  two  successive  orders 
could  be  separated  by  an  arbitrarily  short  time  interval. 
In  this  model,  orders  can  never  be  more  closely  spaced 
than  by  an  interval  of  length  T.   If  T  is  large  enough, 

it  is  possible,  provided  that  there  is  a  sufficiently 
small  range  of  variation  in  the  lead  time,  that  both 
assumptions  hold  simultaneously. 

g.  The  demands  are  Poisson  distributed  with 
mean,  X. 

h.   All  variables  are  treated  as  continuous. 
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3.   System  Two 

This  system  is  identical  to,  but  independent  of, 
System  One. 

D.   ANALYTICAL  SOLUTION  TO  THE  MODEL 

As  was  mentioned  previously,  the  objective  of  the  model 
is  to  minimize  total  system  costs  subject  to  a  constraint 
on  the  number  of  back  orders  per  year. 

The  determination  of  various  annual  costs  of  each 
system  may  be  made  independently  and  placed  into  a  common 
cost  formula.   Instead  of  doing  this  independently,  it  is 
better  to  find  a  main  system  cost  expression,  System  One 
cost  expression  and  System  Two  cost  expression,  and  then 
to  add  them  to  each  other  to  determine  the  total  cost 
formulation  and  then  use  this  formulation  as  the  objective 
function  of  total  system.   The  systems  are  actually  tied 
together  through  the  constraint  on  the  number  of  back 
orders. 

For  the  following  notation,  the  subscript  m  will 
indicate  that  the  variables  being  subscripted  belong  to 
the  Main  Sytem  where  1  and  2  indicate  System  One  and 
System  Two,  respectively. 

The  costs  that  are  of  interest  in  this  model  for  each 
system  are  the  cost  of  placing  an  order,  the  cost  of 
carrying  inventory  and  back  order  costs. 
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1.   Main  System 

As  was  discussed  previously,  the  Main  System  has 
a  transactions  reporting  policy  or  (Q,  r)  model. 

In  a  continuous  review  system,  a  period  is  defined 
as  the  length  of  time  between  the  receipt  of  two  successive 
procurements.   This  time  period  is  a  random  variable. 
Because  the  procurement  lead  times  are  random  variables, 
the  number  of  demands  for  a  fixed  time  are  random  variables; 
and  the  number  of  items  demanded  per  demand  are  also  random 
variables.   The  review  period  is  also  a  random  variable. 
In  the  following  material,  we  itemize  the  costs. 

a.  Procurement  Cost 

A   =  Expected  number  of  demands  per  year 

Q  =  Order  quantity 

A 

pr™-  =  Average  number  of  procurements  per  year 
m 

A  =  Procurement  cost  per  cycle 
m 

A 

tt™-  A  =  Procurement  cost  per  year 

Qm  m 
m 

b.  Inventory  Carrying  Costs 

Because  of  randomness,  it  is  possible  to  accumu- 
late a  large  number  of  back  orders  at  the  end  of  a  cycle. 
To  prevent  this  from  occurring,  it  is  advisable  to  provide 
a  safety  stock  to  buffer  the  system  from  excessive  numbers 
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of  back  orders.   The  safety  stock  is  the  expected  amount 
of  stock  on  hand  when  an  order  arrives.   The  actual  amount 
of  stock  on  hand  when  a  shipment  arrives  is  clearly  random. 

S   =  Mean  value  of  on-hand  stock  when  an  order 
m 

arrives 
After  an  order  arrives,  the  expected  on-hand  inventory  in- 
creases to  Q+S  and  is  reduced  to  a  value  of  S  on  the  average 
just  before  the  next  order  arrives.   Therefore,  the  average 
on-hand  inventory  per  cycle  is: 


(Q+S J    Sm 

m^  m    _m  _  Q  +  c 


2       2    2m 


To  write  S   in  terms  of  the  reorder  point,  r  , 
m  m 

let  us  first  assume  that  the  lead  time  x   is  fixed.   Let 

m 

£  (x;r  )  =  r   -  x  be  the  net  inventory 

at  the  time  an  order  arrives, 

where  x 

is  the  number  of  units  demanded  in  lead  time  t  .   Then, 


S  =  E    [Net  inventory]  =  /   £_(x;r  )  f(x;x  )  dx 

m    t   L  J  J   o      xm    m       m 

m 


where  f(x;x  )  =  Density  function  of  demand  in  time  x_ 
m  J  m 


S   =  /    (r  -x)  f(x;x  )  dx  =  r    /    f(x;x  )  dx 
m   o       m         m        mo         m 

CO 

-  of      f(x''Trm)  dx 
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Sm  =  rm~ym     where  ym  =  Expected  lead  time  demand. 

Then  9?  +  Sm   =  ^   +  r   -  y 
2     m    2     mm 

I  =  Inventory  carrying  charge 

C  =  Cost  of  an  item 

Total  holding  cost/year  =  IC  (r   -  y   +  §?) 

J  m    m    2 

c.   Stockout  Costs 

0  if  x-r  <0 
Let  us  define  n    (x  ;  r  )  =  { 

m  x-r   if  x-r  ^0 

m       m 


where  n    (x;r  )  is  the  number  of  back  orders  per  cycle. 
t      m  c  j. 

m 
If  n  (r  )  ~   expected  number  of  back  orders  per  cycle,  then 

—  00  00 

VW   nxm(X;rm)  h(x)  dx  =  /   ^"V  h(x)  dx 


=    f°°  xh(x)    dx-r      H(r    )    where  H(r    )    =   P  [X>rm] 

r 
m 

and  h(x)  =  marginal  distribution  of  leadtime  demand. 
Therefore,  the  expected  number  of  back  orders/year  = 

f   [  /"  xh(x)  dx  -  rmH(rm)] 
wm   m 

and  the  expected  cost  of  back  orders/year  = 

¥*  [r/"xh(x)  dx-rmH(rm)] 

vm    m 
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All  the  terms  in  the  average  annual  variable 

cost  K   have  now  been  found: 
m 

Km  "   QT  Am   +    IC    hr  +   rm   -   "J    +   Sf    [r'"  *h<x)dx-rm  HtrJ] 

m  *m    m 


2.   System  One  and  System  Two 

It  was  stated  in  the  assumptions  that  System  One  and 
also  System  Two  both  follow  a  periodic  review  policy  with  an 
order  up  to  R  stockage  policy.   Since  System  One  and  System 
Two  are  identical  to  each  other;  the  equations  will  be 
derived  only  for  System  One. 

For  convenience,  a  period  is  assumed  to  be  the  time 
between  the  receipt  of  two  successive  orders  rather  than 
between  the  placement  of  two  successive  orders.   Costs  are 
described  as  follows: 

a.  Ordering  and  Reviewing  Costs 
J  =  Reviewing  cost/cycle 

A  =  Ordering  cost/cycle 

Jl+Al 
Ordering  and  reviewing  cost/year  =  — = 

ll 
where  Tj_  is  the  period  length  defined  in  units  of  years. 

b.  Inventory  Carrying  Cost 

The  expected  net  inventory  just  prior  to  the 
arrival  of  an  order  is  R,  -  y1  -  ^,T  ,  where  u1  =  mean  demand 

during  lead  time. 

The  mean  rate  of  demand  remains  constant  over 
time  and  the  expected  demand  per  period  must  be  the  expected 
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amount  ordered,  i.e.,  ^T,.   If  the  expected  net  inventory 
immediately  after  the  arrival  of  a  procurement  is  R, -y  , 
it  is  therefore  ^  -  y1  -  A  t^  just  prior  to  the  arrival 
of  a  procurement. 

The  expected  unit  years  of  storage  incurred 
per  period  is 

A  T 

Ti  £  (Ri  "  H}   +  \   (Ri  "  »*i  "  W]  ■  Ti[Ri  -  ^1  -  -T1] 

AT 
and  average  inventory  carrying  cost/year  =  IC  [R,  -  y,  -  — = — ] 

c.   Stockout  Costs 

First  we  assume  the  case  where  the  procurement 
lead  time  is  constant  t.   An  order  placed  at  time  t  will 
arrive  in  the  system  at  time  t  +  t ,  and  the  next  procurement 
will  arrive  in  the  system  at  time  t  +  t  +  T,  .   After  the 
order  is  placed  at  time  t,  the  inventory  position  of  the 
system  is  R, .   It  is  necessary  to  compute  the  expected 
number  of  back  orders  occurring  between  t  +  t  and  t  +  T  +  T, . 
A  back  order  will  occur  in  this  period  under  assumption  c   if 
and  only  if  the  demand  in  the  time  period  t  +  T,  exceeds  R. . 
Assumption  c  also  assures  that  after  the  arrival  of  the  order 
placed  at  time  t,  there  will  be  no  remaining  back  orders, 
and  therefore  they  must  all  occur  between  times  t  +  t  and 
t  +  t  +  Tl.   Consequently  the  expected  number  of  back  orders 
incurred  per  period  is 
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/   (x-R-^)  f  (x;  t  +  T1)  dx  where 

f(x;  t  +  T,)  =  Demand  distribution  during  time  t  +  T,  . 

When  lead  time  is  random  with  density  g  (x,) 

with  x  .-   and  x    being  lower  and  upper  limits  respectively 

mm      max     3  rir  c  u 

and  x,  andi-,  the  lead  times  for  the  orders  placed  at  times 
t  and  t  +  T, ,  respectively,  the  expected  number  of  back 
orders  incurred  per  period  is: 

T        T 

max   max 

00 

/  /  R/       (x-^)    f(x;    T2+T1)g(x2)g(T1)dx   dTT2dTl 

min        min  1 


=R/  (x-R^  h(x;  T±)    dx 


^  tti  a  'V 

where  h(x;  T±)    =   /      f (x;  ^2+Tl) g (t2} dT2   which  is  the 

min 

demand  distribution  during  time  x2  +  T,  when  lead  time  is 

a  random  variable  with  density  function  g(xj.   The  average 

number  of  back  orders  incurred  per  year  is; 

E1(R,,.T1)  =  j~-  Rf       (X-R^  h(X;  T^)  dx  and  the  average 

back  order  cost  per  year  equals  to  II  E   (Rp.T^)  and 

n, (r)  =  T  ,  E,  (R, ,T9)  is  the  expected  number  of  back  orders 

per  period. 
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Finally,  the  annual  variable  cost  of  System  One 

is: 

L                   XT 
K-l  =  —  +  IC  [Rx  -  \i1 *j±]    +    J[1   E1  (R1,T1)  and 

likewise  the  annual  variable  cost  of  System  Two  is: 

L  XT 

K2  "  fj  +  IC  fR2  "  ^2  "  4^J  +  n2  E2  (R2'T2) 

where  Li  =  J .  +  A.  . 
1    l 

3.   Objective  Function  and  Constraints 

The  objective  function  of  the  model  consists  of  the 

total  annual  variable  costs  of  each  system.   Therefore; 

K  =  K   +  K,  +  Kn   which  is  equal  to  the  minimization 
m    1    2  n 

of; 

X  q  n  X  L, 

K  =  J£  A m  +   IC  [-5  +  r  -  u   3  +  -^  Hm  (r)  +  =i 

Q    m        2     m    m      Q    m       T, 
wm  m  1 

xt    n  _     l2  x2t21 

+  IC  [Rx  -  p1  -  -^i]  +  ~  nx(r)  +  j±.  +  IC  [R2-y2 2— J 

n2  - 
+  =-  ru  (r) 

x2 

subject  to  ^  nm(r)  +  ~-  n1  (r)  +  ^-  n2  (r)  <b 

ml  2 

where  b  is  the  specified  total  maximum  number  of  back  orders 
per  year  for  the  entire  system. 

The  problem  at  hand  is  to  calculate  the  optimum 
values  of  Qm,  rm,  R-^  T1#  R2 ,  T2-   Since  the  objective 
function  and  the  constraint  are  non-linear  functions  of  the 
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decision  variables,  we  solve  the  problem  using  the  Lagrange 
multiplier  approach. 

4.   Optimum  Values  of  Operating  Variables 

After  including  the  Lagrange  multiplier  in  the 
formulation,  the  new  objective  function  becomes: 

Minimize 

Q  n  x        l, 

T     m  ...__,  r  m  .         -i,   mm-,,,   l 

L  =  tt—  A  +ICI-S-  +  r   -  u   +  —x —  n  (r)  +  =— 

Q-.ni    L  2  m    MmJ     Q    'm   '    T, 

m  ml 

A.T      n   .        L  AT, 

+  ic[Rl  -  Ul  -  _^i]  +  -i  nl(r,  +  2  +  iC[R2  -  u2  -  4^J 

1  " 

+  ji  n2  (r)  -  9[(gS  nm(r)  ♦  i  5  (r)  +  i  n2  (r) )  -  b] 

z  m  1  2 

where  9  is  the  Lagrange  multiplier.   It  is  through  this 
Lagrange  multiplier  that  the  three  systems  are  linked  together 
mathematically . 

The  optimum  values  of  the  unknown  variables  can  be 
found  by  taking  the  derivatives  of  the  objective  function 
with  respect  to  Q  ,  R  ,  R, ,  T,,  R2 ,  T2 ,  0;  equating  them  to 
zero;  solving  the  equations  simultaneously;  and  ensuring  the 
Kuhn-Tucker  conditions  are  satisfied. 

The  derivatives  were  taken  with  respect  to  Q  .  r  , 

E         mm 

R. ,  R2,  6.   A  different  procedure  was  utilized  to  find  T, 
and  T~. 
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The  derivatives: 

T         _    9L                AmAm        IC  Xm     -    ,    ,     ,_ 

Q      ~    30 2~  +   ~2 J-  nm(r)     (n™    "9=0 

m  ^m 

L         -  |i-.   ic   -   J£  H    (r)  (II       -    9)    =    0 

rm        \                    Qm     m  m 
m 


\  =  T£  -  Ic  "  W[  WV    ("i  -  e)  =  o 


■  IC  -  r-  H,    (R„,    T„)     (IL,    -    9)    =    0 


R2         3R2  ~        T2      2    viN2'      2'     l"2 

Le   =  ft  ■  ['q21  Vr)  +  r-  Vr>  +  r-  Vr))  -  »>]  -  o 

vm  1  2 


The  Kuhn-Tucker  conditions: 


L~   ^0       Q  .L_   =  0        Q^O       g  £  0 

Qm  *m  Q  vm  y 

m 


L    *  0       rL=0         r^O       9.g=0 
m  rm 


LR  >,    0       R1'LR  =  °        Rl  *  °      6  <  0 


L-,   ^  0       R0.L_   =  0        R0  ^  0 
R2  2   R2  2 


where  9  =  q~  V<r)  +  t"i  ^l(r)  +  T~~  n2(r)  "  b 


Therefore  the  following  Kuhn-Tucker  conditions  must 
be  satisfied: 
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AmAm    IC    Xm  -  (r)  (II  -  9)  >,    0 

m  m 


X 

IC  -  FT   Hm(r)   (II   -  9)  >,    0 
Q   m      m 
m 


rc  -  —  h1(r1,  tx)  (nx  -  9)  *  0 


IC  -  -i-   H2  (R2,  T2)   (n2  -  9)  >  0 


X  A     IC  Q     A 
Qm       2      Qm  m      m 


X   r 

IC  r   -   mn  m  H  (r)  (II  -  9)  =  0 
m      Q     m       m 
m 


IC  R±   -  ^  k1    {R±,   tx)  (ii1  -  6)  =  0 


IC  R2  -  ^H2  (R2,  T,)   (I2  -  9)  -  0 


Qm 

!> 

0 

r 

m 

> 

0 

Ri 

> 

0 

R2 

> 

0 

n  (r)  +  i  n-,  (r)  +  i  n,(r)  -  b  <  0 


Q   'in 

m 


Tl 


T    2 
X2 


1  - 


1  - 


[^  V<r)  +  ^  Vr)  +^n2(r)  -bi 


=  0 


Q 
S  0 
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From  these  equations  and  conditions  the  optimum 

values  of  the  unknown  variables  can  be  simplified  to  the 

following  form  and  the  exact  values  can  be  found  by 

solving  these  equations  simultaneously: 

IC  Qm 
H  (r)  =       m 


ST  '    (II  -9)  A 

m     m 


On,   =    v/2Xm[ATTi   I    FT      (n      _    9)  ] 
m  m  L    m   "*"      mm _ 

IC 

Hl    (R1'   Tl)    "    (n1-9) 


IC   T 
Hn     (R„,    T0)    =   -7= 


2     ViV2'       2'  (H2-6) 


The  procedure  to  find  T,  and  T2  and  subsequently  R, 
and  Rp  is  through  iteration  and  trial  and  error.   This  is 
not  unreasonable.   In  realistic  cases  other  considerations 
not  modelled  here  usually  dictate  the  length  of  a  period. 
In  most  cases  the  systems  control  a  large  number  of  different 
items  and  the  same  period  length  is  used  for  each  item. 
Thus,  the  period  length  is  usually  some  convenient  calendar 
or  financial  period,  such  as  a  month  or  a  quarter.   The 
procedure  we  use  is  shown  in  the  following  example. 

E.   EXAMPLE  PROBLEM 

1.   Additional  Assumptions 

For  this  example  it  is  assumed  that  each  system  has 
Poisson  arrivals  independent  of  each  other.   The  customers 
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demand  only  one  item  at  a  time.   System  One  and  System  Two 
are  resupplied  only  from  the  main  system  and  the  main  system 
is  resupplied  from  outside  suppliers.   To  be  more  realistic 
it  is  also  assumed  that  the  procurement  lead-times  have  a 
gamma  distribution  with  different  mean  and  variances. 

The  A's  will  be  assumed  daily  demand  or  arrival 
rate  per  day  rather  than  yearly  values. 
2 .   Input  Values 

Before  providing  the  input  values,  some  clarifications 

must  be  provided  regarding  the  demand  at  the  main  system. 

When  it  is  stated  that  Xm  =  4 /day,  A,  =  3/day  and  A2  =  1/day. 

These  represent  the  mean  number  of  demands  which  arrive  each 

day  directly  at  the  respective  systems.   However,  since  all 

demands  at  Systems  One  and  Two  eventually  must  filter  up  to 

the  main  system,  the  cumulative  demand  at  the  main  system 

has  an  expected  value  of  A,  +   \~    +   \      =  8.   The  main  system 

12m 

supports  not  only  Systems  One  and  Two,  but  also  has  its  own 
customer  demands.   The  model  assumes  that  the  demand  at  the 
main  system  is  the   superposition  of  the  direct  demands  at 
the  main  system  and  Systems  One  and  Two.   However,  since  the 
actual  demands  placed  on  the  main  system  in  the  lower  echelon 
are  batched,  the  demand  variability  is  much  greater  than 
would  be  expected  by  the   superposition  process.   We  will 
evaluate  the  seriousness  of  our  assumption  about  the  demand 
process  at  the  main  system  with  the  simulation  model 
described  in  the  next  chapter. 
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are 


In  this  example  the  input  values  for  each  system 


a.  Main  System 

A   =  8/day  (considering  the  individual  customer 
arrivals  and  the  remainder  being  Systems 
One  and  Two  average  demands) 

C   =  $50/item 

I   =  0.23 

m 

Am  =  $500 
m 

b.  System  One 
X±   =    3/day 

C  =  $50/item 
I   =  0.23 

n1  =  $5 

A1   =  $100 
J  =  $100 

c.  System  Two 
\2   =  1/day 

C  =  $50/item 
I   =  0.23 

n2  =  $5 

A2  =  $100 
J2  =  $100 
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The  total  system  is  expected  not  to  exceed 
b  =  175  total  number  of  back  orders  at  the  end  of  a  year. 
3.   Solution  Procedure  and  Results 

It  was  assumed  that  for  each  system,  procurement 
lead  times  are  gamma  distributed  random  variables. 

The  gamma  distribution  is  defined  as: 

t 

/.u\       1     .a-1     8     ... 
g(t)  =  t     e      with 

PL 

Via)    S 

y  =  a*  3  and 

Var  =  a. 6 

We  need  the  lead-time  demand  distribution  for  the 
case  in  which  lead  times  are  gamma  distributed  and  the 
process  generating  demands  is  Poisson.   This  is  derived 
below. 

The  Poisson  distribution  is 

,,,>x   -At 

f (x)  =      ,  e and  the  demand  distribution 

x  J 

during  lead  time  is: 


f(x;x)  =  /   f(x)  g(t)  dt 

o 

.^.x   -At   .a-1   -t/B 
-  r°°    (At)   e     t    e       ,. 
"V   x         T(a)  6a 

A*       rta>l-xe-^+X)tdt 


T(a)  6a  xl  o 
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with 


x  ,1  ,,.a+x.a-l+x    ~^+x^ 

AX  r  (a+x)  rC0(g+A)        t  e 

r(a)    3a  xi       (l+x)a+x0-/  r(a.+  x)  dt 


r   (l+A)a+xta-l+xe-(S+A)t 

0    n^xi dt  - 1    so 


f(^.T\    =  A T  (a+x) 

E  l    '  T;         r(a)    6«   xJ    7T~T7aTx 

(3   +A) 


r  (a+x)     r       A        i  x    r    3       i  ot 
(a)    x!    L  ,1,  „    J        Ll 


(|+A)  f  +  A 


a+x-1  .  4 

(      x      )     [T1-]X    [/-]a 


a+x-1 

=  (      )  (1-P)X  Pa 
x 


a(l-   ^ 


_  a(l-p)  _  "^   1/3  +\'    _  ,  ft 
M l/g Aae 

1/3  +  A 


Var  =  «J±JPL  ,  ^^^   =  Xa3(l+A3) 
p       1/6+X 


In  fact  this  is  a  Negative  Binomial  distribution 


p    1/3+A 
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The  parameters  a  and  8  for  each  system  are: 


am  =  12.8  cc1  =  5.43         a2   =  2.5 

8m  =  3.125  g1  =  4.6  S2  =  6. 


The  main  system  has  a  Negative  Binomial  distributed 
lead  time  demand  with  parameters: 


ym  =  8-  (  12.8)  (3.125)  =  320 
varm  =  8-  (12.S)  (3.125)  [1+8  •  (3 .  125)  ]  =  8320 


System  One  has  the  parameters: 
U1   =  3-  (5.43) •  (4.6)  =  74.934 

var.  =  3-  (5.43)  (4.6)  [1+3(4.6)]  =  1109.0232 


System  Two  has  the  parameters 
U2  =  1- (2.5) -6  =  15 


var2  =  l-(2.5)-6  [1+1(6)]  =  105 


Because  the  negative  binomial  is  computationally 
intractable,  we  use  the  normal  approximation  for  the  calcu- 
lation of  lead  time  demand  probabilities.   The  normal 
distributions  are  assumed  to  have  the  same  mean  and  variance 
as  the  negative  binomial  distributions  they  replace. 

To  solve  this  example  problem,  it  is  necessary  to 
make  an  initial  estimate  for  the  Lagrange  multiplier  9 
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which  must  be  less  than  or  equal  to  zero.   After  the  initial 
estimate  if  the  total  number  of  back  orders  per  year  exceeds 
b,  then  the  absolute  value  of  0  should  be  increased  gradually 
until  the  number  of  back  orders  converges  to  b. 

For  this  problem  6 =  -1.5  works  very  well.   For  the 
main  system  the  formulas  are: 


Qm=  ^VAm  +  Vr)(llnr9)J/lC 


Vr)  =  ( Vmm  =  P  tx^J 


Q   is  calculated  to  be  550  and 


H 
m 


(r)  =  P  ^*J    =  5(g+{°M32920  =  °-3332455216 
From  the  inverse  standard  normal  distribution  TI-59 
calculator  program  (Appendix  G) ,  the  reorder  level  at  the 
main  system   is  found  to  be: 


r   =  a  .(0.4305333395)  +  y   =  359.2706827 

mm  Hm 


r  -y 

n  (r)  =  (y  -r  )  H  (r)  +  a   0  (— — -) 
'm       Hm  m   m       ma   ; 

m 


where  0   is  the  functional  value  of  standard  normal  distribu- 
tion at  0.4305333395 

n  (r)  =  20.08139334.   Using  this  value  in  the  Qm 
m  m 


formula  yields 
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Qm  =  /2v2gg(j  [500  +  20.08139334  (5  +  TTBT    -co    OA 
(0.23). 50 =  568.84 

which  is  not  equal  to  the  first  estimated  value.   If  Q   is 

m 

taken  56  9  then: 


r   =  a       (0.399071201)  +  u      «  356.4008941 

mm  m 


fL(r)    =     21.05458646 
m 

Using  this   value    in   the   Q      formula  yields 

m         -* 


Q   =  /2.2920  [500  +  21.05438646  (5  +  1.5)'=  568.69 

(0.23)-50 

which  is  very  close  to  the  initial  Q   value.   A  summary  of 

the  results  for  the  main  system  follows: 

Qm        =  569 
m 

S     =  36.400891 

m 

n  (r)  = 21.05438646/period 

rL(r)  = 108.047115/year 
m  J 

K    =  $6796.50/year 
m 

To  solve  for  the  optimum  values  for  Systems  One  and 

Two  a  different  procedure  is  followed.   The  total  annual  cost 

of  Systems  One  and  Two  is  a  convex  function  of  the  period  length 

T.   For  different  values  of  T  there  are  different  total  annual 


33 


cost  values.   The  minimum  of  these  values  is  the  optimum 
total  annual  variable  cost  and  the  corresponding  T  and  R 
values  are  the  optimum  operating  values.   A  TI-59  program 
was  written  to  perform  the  line  search  for  the  best  value 
of  T.   The  program  is  found  in  Appendix  G.   User  information 
and  the  features  of  the  program  are  also  in  the  same  appendix, 
This  program  evaluates  the  R  value,  safety  stock,  back 
orders  per  period,  back-orders  per  year,  annual  reviewing 
and  ordering  cost,  annual  holding   cost,  annual  back  order 
cost,  and  finally,  the  total  annual  cost. 

Using  this  program  the  optimum  values  for  Systems 
One  and  Two  are: 


Rx      =  306.809389 


T,      =  2.39  months  (72.6958  days) 


S1      =  13.78788901 


rutr)      =    8.670289086/period 


n1(r)      =  43.53283223 


K,      =  $2634.41/year 


R2  =  134.9807003 

T2  =  4.07  months  (123.7958333  days) 

s2  =  0 

n2 (2)  =  8.13480099/period 
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n2(r)   =  23.98467122/year 
K2      =  $1377.56/year 

The  total  back  orders  per  year  are: 

nm(r)  +  nx(r)  +  n2(r)  =  108.047  +  43.532  +  23.984 

=  175.56 

which  is  the  same  as  b  =  175.56. 

The  total  annual  variable  cost  is: 

K  =  K   +  K,  +  K0  =  6796.50  +  2634.41  +  1377.56 
m    1    2 

K  =  $10808.47. 
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III.   COMPUTER  SIMULATION  FOR  MULTT -ECHELON  MULT I -LOCATION 
SINGLE  ITEM  INVENTORY  SYSTEM 

To  check  the  analytical  results,  a  computer  simulation 
was  written  which  uses  the  same  operating  assumptions  and 
input  parameters  that  were  made  for  the  analytical  model. 
This  model,  however,  simulates  the  real  world  more  accu- 
rately since  some  of  the  simplifying  assumptions  required 
to  obtain  analytical  results  were  not  necessary  in  the 
simulation.   Also,  the  demands  placed  at  the  main  system 
from  the  lower  echelon  systems  were  batched  as  in  the  real 
world. 

A.   DESCRIPTION  OF  THE  SIMULATION  MODEL 

As  in  the  analytical  case,  there  are  three  systems  in 
the  simulation  model;  the  main  system,  System  One  and  System 
Two.   The  flow  charts  of  this  program  are  in  Appendix  A. 

1.   Main  System 

This  system  uses  a  continuous  review  policy.   When 
the  stock  level  reaches  the  reorder  point,  it  orders  the 
amount  Q.   The  decision  variable  for  reordering  is  the  in- 
ventory position.   When  an  order  is  placed,  the  inventory 
position  increases  by  an  amount  of  Q.   If  the  new  inventory 
position  is  less  than  the  reorder  point,  the  system  places 
an  additional  order  for  another  amount  of  size  Q.   Then  the 
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inventory  position  increases  with  one  more  Q.   This  con- 
tinues until  the  inventory  position  exceeds  r.   The  order 
policy  is  thus  (nQ,r).   N  is  the  smallest  integer  that  will 
make  the  inventory  position  higher  than  the  reorder  point. 

There  are  independent  customer  arrivals  to  the 
main  system  and  also  group  demands  that  are  placed  by 
System  One  and  System  Two  when  those  systems  replenished 
their  stocks  at  the  end  of  their  periods. 

The  number  of  units  demanded  per  requisition  is 
one.   The  number  of  units  demanded  for  resupply  to  the  lower 
echelon  systems  is  random  depending  on  the  demands  at  the 
lower  echelon  systems  and  the  parameters  R, ,  R_ .   The 
program  is  capable  of  allowing  geometrically  distributed 
quantities  demanded  per  requisition.   Back  orders  at  the  main 
system  are  satisfied  by  filling  the  back  orders  to  individual 
customers  first  followed  by  filling  any  back  orders  due  to 
System  One  and  System  Two  on  a  first-come,  first-served 
basis. 

When  a  demand  occurs  from  the  lower  echelon  systems 
for  which  there  is  not  sufficient  on-hand  inventory  at  the 
main  system,  the  maximum  amount  is  filled  and  the  rest  is 
put  into  the  back  order  queue.   As  soon  as  a  shipment 
arrives  at  the  main  system,  the  back  orders  are  filled. 
The  times  between  the  customer  arrivals  are  independent 
of  each  other  and  exponentially  distributed.   The  lead 
times  are  also  independent  of  each  other  and  gamma  distributed. 
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2.   System  One  and  System  Two 

Both  systems  operate  identically.   The  only  differ- 
ence might  be  in  the  values  of  the  system  parameters.   They 
have  the  same  periodic  review  policy  which  is:   order  up  to 
R  at  each  review  time.   The  program  is  also  capable  of  making 
a  decision  at  each  review  time  regarding  the  placement  of 
an  order  if  the  inventory  position  is  less  than  or  equal  to 
a  threshold  value,  is  an  (r,R)  policy.   If  r  is  taken  to  be 
R-l,  then  the  system  orders  up  to  R  at  each  review  time  even 
if  there  is  only  one  demand  in  a  period. 

The  times  between  customer  arrivals  are  independent 
of  each  other  and  exponentially  distributed.   The  number  of 
units  demanded  per  requisition  is  one.   The  lead  times  are 
independent  of  each  other  and  gamma  distributed. 

B.   HOW  TO  USE  THE  PROGRAM 

Two  simulations  are  given  in  Appendix  C  and  Appendix  D. 
The  programs  differ  primarily  in  the  type  of  output  that 
is  generated.   The  program  given  in  Appendix  C  produces  a 
Versatec  plot  output  showing  the  inventory  position  for 
the  main  system  and  the  lower  echelon  systems.   A  simulated 
period  of  up  to  4  or  5  years  can  be  run  with  a  Class  K  job. 

On  the  Versatec  output  for  each  system,  there  are  two 
plots:   one  plot  has  a  small  triangle  at  each  point  repre- 
senting the  inventory  position.   The  second  plot  represents 
the  net  inventory.   The  first  figure  shows  the  main  system 
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values,  the  second  figure  shows  the  System  One  values  and 
the.  third  figure  shows  System  Two  values  all  plotted  versus 
time. 

The  second  program  can  be  used  to  simulate  operations 
over  arbitrarily  long  periods  of  time.   There  are  no  dimen- 
sion restrictions.   This  program  gives  the  net  inventory 
at  the  end  of  each  period  for  all  the  systems;  total  demand 
in  the  periods  of  both  System  One  and  System  Two;  demand 
during  a  lead  time  for  the  main  system;  average  on-hand 
inventory;  and  average  number  of  items  short  per  day  for 
each  system. 

The  input  required  by  the  programs  is  described  by  the 
variable  definition  list  given  in  Appendix  H.   The  programs 
use  the  IMSL  subroutines  for  random  number  generation. 

C .   RUNS 

1.   Starting  Conditions 

There  were  two  different  starting  conditions  entered 
for  the  main  system.   In  the  first  case  the  inventory  posi- 
tion and  net  inventory  were  set  initially  to  equal  the  order 
quantity.   In  the  second  case,  the  inventory  position  and 
the  net  inventory  were  set  equal  to  the  reorder  point  plus 
the  order  quantity.   For  Systems  One  and  Two  the  inventory 
positions  and  net  inventories  were  set  initially  at  R. . 
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2.   Results  of  Runs 

In  all  runs  the  length  of  time  was  10  years.   The 
average  yearly  results  were  obtained  by  dividing  the  results 
for  the  10  year  period  by  10. 

Using  the  policy  parameters  determined  by  the  mathe- 
matical model  and  the  first  set  of  starting  conditions  the 
following  simulated  results  were  obtained. 

a.  Main  System 

The  number  of  orders  in  10  years  =  51. 

The  number  of  back  orders  in  10  years  =  3026. 

The  average  on-hand  inventory  over  10  years  = 

313.026. 

The  average  safety  stock  when  an  order  arrives  = 

56.54  units. 

The  average  annual  costs  are: 

Ordering  costs:      $2,550.00 

Holding  costs:        3,599.80 

Back  Order  Costs:     1,513.00 

Total  Cost    $7,662.80 

b.  System  One 

Number  of  back  orders  in  10  years  =  587. 

The  average  on-hand  inventory  over  10  years  = 

114.482  units. 
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The  average  annual  costs  are: 

Ordering  and  Reviewing  Costs:  $1,004.18 

Holding  Costs:  1,316.54 

Back  Orders  Costs:  2  93.50 


Total  Costs  $2,614.22 

c.   System  Two 

Number  of  back  orders  in  10  years  =  287. 

The  average  on-hand  inventory  over  10  years  = 

57.618  units. 

The  average  annual  costs  are: 

Ordering  and  Reviewing  Costs:   $   589.68 

Holding  Costs:  662.61 

Back  Orders  Cost:  143.50 

Total  Costs  $1,395.79 

(1)   Total  results 

The  total  number  of  back  orders  for  the 

whole  system  per  year  is  302.6  +  58.7  +  28.7  =  390. 

The  total  annual  variable  cost  for  the 

whole  system  is  $7,662.80  +  $2,614.22  +  $1,395.79  =  $11,672.81. 

The  model  was  run  again  with  the  second  set 

of  starting  conditions  and  the  same  set  of  policy  parameter 

values.   The  simulation  results  are  summarized  below: 

(a)   Main  System 

The  number  of  orders  in  10  years  =  50. 

The  number  of  back  orders  in  10  years  = 

1665. 
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The  average  on-hand  inventory  over 

10  years  =  346.937. 

The  average  safety  stock  when  an  order 

arrives  =  77.46. 

The  average  annual  costs  are; 

Ordering  Costs:       $2,500.00 

Holding  Costs:         3,989.78 

Back  Order  Costs:        832.50 

Total  Costs:      $7,322.28 

(b)  System  One 

Number  of  back  orders  in  10  years  =  582. 

The  average  on-hand  inventory  over 

10  years  =  118.696  units. 

The  average  annual  costs  are; 

Ordering  and  Reviewing  Costs:  $1,004.18 

Holding  Costs:  1,365.00 

Back  Orders  Costs:  291.00 

Total  Costs:  .  $2,660.18 

(c)  System  Two 

Number  of  back  orders  in  10  years:   314. 
The  average  on-hand  inventory  over  10  years 
57.401  units. 
The  average  annual  costs  are; 
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Ordering  and  Reviewing  Costs:  $  589.6  8 

Holding  Costs:  660.11 

Back  Orders  Costs:  157.00 

Total  Costs:  $1,406.79 

(2)   Total  Results 

The  total  number  of  back  orders  for  whole  system 

per  year  is  166.5  +  58.2  +  31.4  =  256. 

The  total  annual  variable  cost  for  whole  system 

is  $7,322.28  +  $2,660.18  +  $1,406.79  =  $11,389.25. 

3.   Comparison  of  Both  Starting  Conditions 

It  is  obvious  that  if  both  results  are  compared  there  is 

a  considerable  decrement  on  annual  back  orders  but  not  a  great 

difference  in  annual  variable  costs. 

390-256   _  n  -.-.. 
390    ~  °-3436 

11672.81  -  11389.25 


11672.81 


=  0.02429 
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IV.   COMPARISON  OF  ANALYTICAL  AND  COMPUTER  SIMULATION  RESULTS 

By  comparing  the  analytical  results  with  the  simulation 
results,  it  is  possible  to  evaluate  the  analytical  model. 
If  the  measures  of  effectiveness  predicted  by  the  analytical 
model  are  reasonably  close  to  those  generated  by  the  simula- 
tion, there  is  support  for  the  analytical  results.   In  the 
table  below  we  compare  the  main  measures  of  effectiveness: 
back  orders  and  costs. 


Analytical  result 

Simulation  result 
(For  both  starting 
conditions) 


Total  Back  Orders   Total  Variable 
Per  Year        Costs  Per  Year 


175.56 
390.00 
256.1 


$10,808.17 
$11,672.81  or 
$11,389.25 


This  table  shows  that  the  average  yearly  back  orders  are 
considerably  higher  than  what  is  estimated  by  the  mathematical 
solution.   In  fact,  the  total  number  of  back  orders  generated 
in  the  simulation  does  not  even  meet  the  constraint.   Since 
the  solutions  produced  by  the  simulation  model,  especially 
with  respect  to  the  number  of  back  orders  do  not  agree  with 
the  mathematical  model,  this  suggests  that  some  of  the 
assumptions  made  in  deriving  the  mathematical  results  are  not 
reasonable.   The  major  differences  between  the  simulation 
results  and  the  analytical  results  are  found  in  the  measures 
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of  effectiveness  at  the  main  system.   The  results  for  the 
lower  echelon  track  reasonably  closely. 

The  costs  are  also  a  little  higher  than  what  is 
expected  but  much  closer  percentage  wise. 

As  an  explanation  for  the  large  differences  observed 
in  the  number  of  back  orders  at  the  main  system  with  the 
simulation  and  analytical  models,  let  us  reexamine  the 
assumptions  we  made  that  affect  back  orders.   Back  orders 
results  from  inadequate  amounts  of  safety  stock  to  protect 
the  inventory  system  against  excessively  large  numbers  of 
demands  in  a  lead  time.   The  safety  stock  is  manipulated 
by  control  of  the  reorder  point  r.   After  the  reorder  point 
is  hit  and  an  order  placed,  the  system  is  totally  at  the 
mercy  of  the  demands  that  occur  during  the  lead  time.   If 
the  variance  of  lead  time  demand  is  underestimated,  then 
large  stockouts  will  occur,  even  if  the  mean  lead  time 
demand  is  estimated  accurately.   The  lead  time  demand  is 
affected  by  two  random  quantities:   (1)   the  demand  distri- 
bution and  (2)  the  lead  time  distribution. 

For  purposes  of  making  the  mathematics  tractable,  we 
assumed  in  our  analytical  model  that  demands  at  the  main 
system  flowed  in  at  a  smooth  continuous  rate  A  which  was 
taken  to  be  the  sum  of  the  rates  of  demand  incurred  directly 

at  the  main  system  A   and  those  which  occurred  at  the  lower 

■*      m 

echelon  systems  A,  and  X-.   The  lead  times  were  assumed  to 
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be  normally  distributed  with  mean  equal  to   At  and  variance 
\(1-Ag)  where  x  is  the  mean  value  of  the  lead  time. 

Because  of  the  relatively  high  demand  rates  used  in  the 
example  runs,  the  normal  assumption  should  be  justified  by 
the  Central  Limit  Theorem.   However,  as  is  illustrated  by 
the  Versatec  plot  output  (See  Appendix  F)  the  demands  at  the 
main  system  are  far  from  smooth.   What  happens  is  the  demands 
which  arrive  directly  at  the  main  system  cause  the  inventory 
position  to  drop  off  smoothly.   However,  when  the  replenish- 
ment orders  from  the  lower  echelon  are  received,  large  drops 
occur  in  the  inventory  position  of  the  main  system.   Recall 
that  the  lower  echelon  systems  order  in  batches  once  each 
period  from  the  main  system.   If  the  main  system  simply  tries 
to  average  out  demands  (as  assumed  by  the  analytical  model) , 
it  will  sometimes  have  very  large  amounts  of  excess  stocks 
when  shipments  arrive  and  sometimes  very  large  numbers  of 
back  orders.   The  high  variance  in  lead  time  demand  caused 
by  the  irregular  demand  actually  seen  at  the  main  system 
causes  the  problem.   Since  all  demands  eventually  flow 
through  the  main  system  and  this  is  assumed  by  our  analytical 
model,  the  problem  cannot  lie  in  the  value  used  for  the 
mean  lead  time  demand. 

Let  us  explore  further  what  happens  when  reorders  are 
triggered  at  the  main  system.   Demands  directly  at  the  main 
system  eat  away  smoothly  at  the  inventory  position.   Then  a 
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very  large  quantity,  say  X,  is  demanded  by  one  of  the  lower 
echelon  systems.   There  is  a  very  good  chance  that  the 
large  order  placed  by  the  lower  echelon  system  will  trigger 
a  reorder  by  the  main  system.   However,  if  the  demand  causes 
a  large  overshoot  of  the  reorder  point,  the  main  system  may 
have  much  less  stock  to  live  off  of  until  the  shipment 
arrives.   For  example,  suppose  the  inventory  position  is 
IP  =  347,  the  reorder  point"  is  r  =  300,  and  System  One  places 
it's  resupply  order  for  200  units.   An  order  will  be  placed 
by  the  main  system  but  instead  of  having  300  units  to  keep 
it  going  until  the  order  is  received,  it  will  have  only 
147  units.   It  is  clear  that  if  300  is  the  amount  of  stock 
needed  to  provide  reasonable  protection  against  demands  in 
a  lead  time,  large  numbers  of  back  orders  would  be  expected. 
Effectively,  in  the  example  above,  the  reorder  level  was  not 
r  =  300  but  r  =  147  and  the  safety  stock  negative.   The 
impact  of  this  surge  in  demand  caused  by  the  batching  of 
demands  received  from  the  lower  echelon  systems  is  to  reduce 
the  "effective"  reorder  point  from  the  value  r  to  a  value 
r'<r.   Because  the  actual  demand  distribution  witnessed  b, 
the  main  system  is  difficult  to  describe  mathematically, 
the  actual  value  of  r'  cannot  be  determined  analytically. 
However,  it  is  clearly  less  than  r  and  may  be  much  more  so. 
In  the  next  chapter  we  describe  a  policy  modification 
for  operation  of  the  main  system  that  was  suggested  by  the 
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observations  above.  The  modification  attempts  to  allow 
the  main  system  to  anticipate  the  surge  of  demands  that 
will  be  received  by  the  lower  echelon  systems. 
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V.   AN  ALTERNATIVE  SOLUTION 

A.   DESCRIPTION  OF  ALTERNATIVE  PROCEDURE 

Due  to  the  reasons  mentioned  in  Chapter  IV,  the  number 
of  back  orders  found  by  simulation  were  much  higher  than 
the  number  of  back  orders  predicted  by  the  mathematical 
model.   The  question  is  how  to  run  the  multi-echelon  system 
so  that  the  large  number  of  back  orders  seen  earlier  can 
be  reduced. 

Obviously,  what  needs  to  be  done  is  to  reduce  the  impact 
of  the  very  large  demands  that  occur  when  the  lower  echelon 
systems  place  their  resupply  orders.   With  modern  day 
communication  and  data  systems,  it  would  be  feasible  to 
allow  the  main  system  to  "see"  every  demand  that  occurs 
anywhere  in  the  system.   If  the  main  system  is  given  the 
visibility,  it  will  be  able  to  take  action  to  get  the  stock 
on   its  shelves  in  anticipation  of  the  large  demands  from 
the  lower  echelon  systems.   The  main  system  can  do  this  if 
it  uses  as  its  reorder  point  the  pseudo  inventory  position 
which  is  like  the  inventory  position  except  that  it  decreases 
only  when  direct  customer  demands  are  encountered  at  any 
of  the  three  systems.   The  pseudo  inventory  position  is  un- 
affected by  the  batch  replenishment  demands  placed  by  the 
lower  echelon  systems.   For  example,  if  a  customer  requests 
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a  unit  from  System  One,  the  pseudo  inventory  position  at  the 
main  system  decreases  by  one.  However,  if  System  One  places 
an  order  for  200  units  at  the  main  system,  the  pseudo  inven- 
tory position  does  not  change. 

This  new  policy  is  referred  to  in  this  thesis  as  the 
"early  warning  policy."   Clearly,  the  pseudo  inventory  posi- 
tion will  always  reach  the  reorder  point  before  the  inventory 
position.   Therefore,  orders  will  always  be  placed  earlier 
and  consequently,  the  number  of  back  orders  should  decrease. 
The  price  paid  will  be  in  terms  of  extra  holding  costs.   The 
results  of  a  simulation  using  the  early  warning  policy  should 
be  more  nearly  like  those  of  the  mathematical  model  since, 
in  effect,  the  mathematical  model  makes  the  assumption  of 
early  warning.   Implicitly,  assumption  of  demand  which  is 
the  superposition  of  the  direct  demands  occurring  at  the 
main  system  and  Systems  One  and  Two  is  equivalent  to  the 
"early  warning  assumptions." 

In  the  next  section,  results  are  given  of  a  simulation 
of  the  multi-echelon  system  with  early  warning.  The  flow- 
chart of  the  simulation  model  is  given  in  Appendix  B.  The 
actual  FORTRAN  computer  program  is  given  in  Appendix  E. 

B.   COMPUTER  SIMULATION  RESULTS  FOR  EARLY  WARNING  POLICY 

The  same  starting  conditions  and  input  parameters  used 
in  the  previous  simulation  run  were  utilized  here.   The 
results  are  summarized  as  follows: 
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1.  Main  System  Results 

The  number  of  orders  in  10  years  =  50. 

Total  number  of  back  orders  in  10  years  =  126. 

The  average  on-hand  inventory  =  4  88.093. 

Average  safety  stock  =  210.96. 

The  costs  are: 

Ordering  costs  =  $2,500.00 

Holding  costs   =   5,613.00 

Back  order  costs     63.00 

Total  costs   $8,176.07 

2.  System  One  Results 

Total  number  of  back  orders  is  10  years  =  665. 

The  average  on-hand  inventory  over  10  years  =  117.846 

The  costs  are: 

Ordering  and  reviewing  costs  =  $1,006.18 

Holding  Costs  =   1,355.23 

Backorder  costs  =     332. 50 

Total  costs  $2,691. 91 

3.  System  Two  Results 

Total  number  of  back  orders  in  10  years  =  312. 

The  average  on-hand  inventory  over  10  years  =  58.552. 

The  costs  are: 

Ordering  and  reviewing  costs  =  $   589.68 

Holding  costs  =  673.35 

Back  orders  costs  =  156. 00 

Total  costs  $1,419.03 

C.   COMPARISON  OF  EARLY  WARNING  POLICY  SIMULATION  RESULTS 
WITH  ANALYTICAL  AND  FIRST  SIMULATION  RESULTS 

The  comparison  will  be  done  with  respect  to  total 
yearly  back  orders  and  total  variable  system  costs  (con- 
sidering all  three  system  entities) . 
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Entire  System 


Number  of  Back    Total 
Orders  Per  Year    Cost 


Analytical  Result  175.56     ,  $10,808.47 

First  Simulation  Results       256.1        11,389.25 

Early  Warning  Simulation 
Results  110.3         12,287.01 

Since  the  effects  of  implementing  the  early  warning  policy 
will  be  observed  primarily  at  the  main  system,  we  also 
produce  the  results  obtained  for  the  main  system  individually 

Main  System 


Numbers  of  Back    Total 
Orders  Per  Year     Cost 


Analytical  Result  108.05       $6,796.50 

First  Simulation  Result        166.50        7,322.28 

Early  Warning  Simulation 
Results  12.6  8,176.07 

This  table  shows  the  differences  better  than  the  first  one. 
The  number  of  back  orders  decreases  88.3  percent,  simultaneously 
as  the  costs  increase  about  10.44  percent.   The  reason  for  the 
higher  cost  is  because  the  increment  in  safety  stock  increases 
the  carrying  costs. 
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VI.   CONCLUSIONS 

The  differences  between  the  results  of  the  mathematical 
and  the  simulation  model  can  be  explained  largely  as  a  result 
of  the  assumptions  made  about  the  demand  process.   In  the 
main  system  it  was  assumed  that  the  demand  was  smooth  (the 
superposition  of  three  Poisson  processes) ,  but  in  the  simu- 
lation model  the  actual  demand  at  the  main  system  was  as  it 
would  be  in  actual  practice.   There  were  batches  of  demand 
placed  by  System  One  and  System  Two  in  addition  to  the  indi- 
vidual customer  demands  directly  at  the  main  system.   These 
demand  batches  in  fact  increased  the  variance  in  lead  time 
demand  beyond  that  modelled.   This  explains  why  many  more 
back  orders  were  generated  in  the  simulation  model  than 
what  was  predicted  by  the  mathematical  model. 

The  simulation  model  developed  in  this  thesis  is  useful 
for  making  comparisons  and  examining  the  effects  of  policy 
changes  or  parameter  changes.   Moreover,  it  is  one  of  the 
best  ways  to  check  the  reasonableness  of  all  of  the  simpli- 
fying assumptions  made  in  order  to  obtain  analytical  solutions 

The  mathematical  results  described  in  this  thesis  do 
not  adequately  determine  the  reorder  point.   The  predicted 
number  of  stockouts  is  much  less  than  the  simulated  numbers. 
As  explained  earlier,  this  is  probably  due  to  the  assumptions 
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made  in  the  model  about  the  variance  of  lead  time  demand. 
Further  study  is  needed  to  determine  what  could  be  done 
in  the  analytical  model  to  better  approximate  what  happens 
in  actual  practice. 

The  early  warning  policy  discussed  in  this  thesis  did 
provide  for  great  reductions  in  the  number  of  stockouts 
system  wide.  Since  stockouts  are  probably  the  most  important 
consideration  in  military  supply  systems,  the  early  warning 
policy  is  recommended,  even  though  the  holding  costs  are 
larger.   Additional  study  is  required  to  see  if  the  reorder 
level  in  the  early  warning  policy  can  be  reduced  substantially 
from  the  value  determined  by  the  mathematical  model.   Pre- 
liminary evidence  is  that  the  reorder  level  can  be  reduced 
significantly  (25  percent  or  so)  without  generating  excessively 
many  back  orders  if  the  early  warning  policy  is  used. 

Our  simulation  models  allow  us  to  view  the  effect  of 
changes  but  they  cannot  be  used  to  optimize  the  values  of 
the  policy  parameters.   For  that  objective,  additional  work 
in  the  mathematical  modelling  area  is  required. 

In  this  thesis  we  have  tried  to  model  a  multi-echelon 
inventory  system  analytically,  by  linking  together  the 
individual  echelons  and  locations  through  a  single  objective 
function  and  a  constraint  on  backorders  system  wide.   We, 
knowingly,  were  making  various  simplifying  assumptions  to 
facilitate  the  derivation  of  solutions.   As  reported  above, 
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the  resulting  solution  for  the  reorder  level  at  the  main 
system  led  to  many  more  back,  orders  than  predicted.   The 
other  solutions;  R  ,  R2 ,  Q,  T. ,  and  T_  appear  to  be  satis- 
factory. 

In  order  to  model  adequately  the  multi-echelon  system, 
it  will  probably  be  necessary  to  build  into  the  determina- 
tion of  the  reorder  lead  at  the  main  system  the  values  of 
the  parameters  R,  and  R~  at  the  lower  echelon  systems. 
This  will  be  the  only  way  to  accurately  describe  the  actual 
demand  process  that  is  observed  at  the  main  system.   We 
recommend  future  work  in  this  area. 
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APPENDIX  A 
FLOWCHARTS  OF  FIRST  SIMULATION  MODEL 

Subroutine  One  =  Ship  arrivals  to  System  One 
Subroutine  Two  =  Ship  arrivals  to  System  Two 
Subroutine  Three  =  Ship  arrivals  to  Main  System 
Subroutine  Four  =  Periodic  review  of  System  One 
Subroutine  Five  =  Periodic  review  of  System  Two 
Subroutine  Six  =  Shipment  arrival  to  System  One 
Subroutine  Seven  =  Shipment  arrival  to  System  Two 
Subroutine  Eight  =  Shipment  arrival  to  Main  System 
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APPENDIX  B 
Flowcharts  of  Early  Warning  Simulation  Model 


Subroutine  One  =  Ship  arrivals  to  System  One 
Subroutine  Two  =  Ship  arrivals  to  System  Two 
Subroutine  Three  =  Ship  arrivals  to  the  Main  System 
Subroutine  Four  =  Periodic  review  of  System  One 
Subroutine  Five  =  Periodic  review  of  System  Two 
Subroutine  Six  =  Shipment  arrival  to  System  One 
Subroutine  Seven  =  Shipment  arrival  to  System  Two 
Subroutine  Eight  =  Shipment  arrival  to  the  Main  System 
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APPENDIX    C 
Simulation   Program;    Versatec   Plotter 
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APPENDIX    D 
Simulation   Program;    No  Versatec   Plotter  Output 
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APPENDIX   E 
Early  Warning   Simulation   Program 
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APPENDIX  F 
Versatec  Output  of  the  Program  in  Appendix  C 
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APPENDIX  G 


TI-59  Calculator  Program  for  Analytical  Solutions  for 
Periodic  Review  Systems 


This  program  was  written  for  probabilistic  periodic 
review  inventory  models  having  gamma  distributed  lead  times 
and  Poisson  arrivals  and  uses.   Normal  distribution  for 
demand  during  lead  time  instead  of  Negative  Binomial. 
Input  Requirements 

The  following  variables  should  be  stored  in  the  regis- 
ters shown  before  the  variables. 

STO  01  =  A  (arrival  rate  per  day) 
STO  02  =  C 
STO  03  =  I 
STO  04  =  n 
STO  05  =  J 
STO  06  =  A 
STO  07  =  a-1 
STO  08  =  S 
This  program  also  requires  the  TI-59  applied  statistics 
module . 
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APPENDIX  H 
Variable  Definitions  for  Simulation  Programs 

IQ  =  Net  inventory 

IP  =  Inventory  position 

IPP  =  Second  inventory  position  for  early  warning  system 

IR  =  Order  quantity 

IS  =  Reorder  level 

ID  =  Number  of  items  demanded  by  systems 

IF  =  Number  of  items  demanded  per  demand  by  ships 

QU  =  Ship  queue 

QS     =  Lower  echelon's  group  demand  queue  in  the  main 
system 

IG     =  Amount  of  demand  for  each  demand  waiting  in  the 

queue 

IA     =  Amount  of  shipment  arrived 

IM     =  Index.   If  it  is  equal  to  1,  that  means  that  the 
demand  waiting  in  the  main  system  queue  to  be 
filled  belongs  to  System  One. 

TW  =  Total  waiting  time 

TOH  =  Total  average  on-hand  inventory 

T  =  Length  of  a  period 

EVENT  =  This  indicates  the  subroutines 

IEVENT  =  This  indicates  time  of  subroutines  scheduled 

X  =  Net  inventory  variable  for  Versatec  plotter 

V  =  Inventory  position  variable  for  Versatec  plotter 

Y  =  Time  for  Versatec  plotter 
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WK  =  Work  space  for  geometric  random  variable 

WS  =  Work  space  for  geometric  random  variable 

WZ  =  Work  space  for  geometric  random  variable 

S  =  Exponential  random  number 

IK  =  Geometric  random  number 

WW  =  Increment 

I      =  Indicates  the  number  of  ships  waiting  in  the  ship 
queue  at  System  One 

J      =  Indicates  the  number  of  ships  waiting  in  the  ship 
queue  at  System  Two 

L      =  Indicates  the  number  of  ships  waiting  in  the  ship 
queue  at  the  main  system 

K      =  Indicates  the  number  of  demand  batches  waiting 
in  the  group  demand  queue  at  the  main  system 

SS     =  Time  indicator 

IQQ    =  Counter  for  ship  arrivals  to  indicate  the  number 
of  items  demanded  per  period 

ML     =  Multiplier  for  the  number  of  batches  of  demand  to 
be  ordered  from  outside  supplier 

KA     =  Indicates  the  last  change  on  net  inventory  of  the 
main  system  in  order  to  get  the  number  of  items 
demanded  in  a  lead  time 

AK  =  Time  of  last  change 

A  =  Scale  parameter  for  a  gamma  distribution 

B  =  Shape  parameter  for  a  gamma  distribution 

XM  =  Ship  arrival  rate  per  day 

P  =  Probability  of  success  for  geometric  distribution 
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